Superradiant Solid in Cavity QED Coupled to a Lattice of Rydberg Gas 
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We study an optical cavity coupled to a lattice of Rydberg atoms, which can be represented by a 
generalized Dicke model. We show that the competition between the atomic interaction and atom- 
light coupling induces a rich phase diagram. A novel "superradiant solid" (SRS) phase is found, 
where both the superradiance and crystalline orders coexist. Different from the normal second order 
superradiance (SR) transition, here both the Solid-1/2 and SRS to SR phase transitions are first 
order. These results are confirmed by the large scale quantum Monte Carlo simulations. 
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Introduction - The study of quantum many-body prob- 
lems and quantum phase transitions (QPT) has attracted 
great interest and is currently one of the main issues in 
the condensed matter community [1]. In the past decade, 
the successful control of the interaction strength and di- 
mensionality of ultracold quantum gases has made it pos- 
sible to explore many interesting physical phenomena [2] . 
For example, the observation of the superfluid-Mott in- 
sulator transition [3], a Tonks gas [4], the BEC to BCS 
crossover [5], and the Kosterlitz-Thouless phase transi- 
tion [6] have all been established in ultra-cold quantum 
gases. 

More recently ultracold atoms have been combined 
with cavity QED [7, 8] to study atom-light coupled many- 
body problems, which give rise to new phenomena. In 
particular, when the two-level atom gas is coupled to a 
cavity [9] the coherent non-local atom-light interaction 
supports the famous superradiance (SR) phase transi- 
tion in the Dicke model (DM) [10, 11]. This SR phase is 
formed by the condensation of atom-light coupled polari- 
tons [12] and breaks the U(l) symmetry. However, the 
interactions between atoms, which may bring new phe- 
nomena or induce novel QPTs, are not considered in the 
DM. 

Atomic interactions promise to give interesting new ef- 
fects and can be implemented via Rydberg atoms. The 
unique properties of strong dipole-dipolc interactions and 
quite long lifetimes have made them a powerful tool 
for the implementation of coherent blockade effects and 
quantum information [13]. Especially, the successful 
trapping of Rydberg atoms in a ID optical lattice [14] 
has stimulated the study of many-body quantum systems 
such as spin system [15, 16] and dynamical crystallization 
or melting of ultracold atoms [17, 18]. 

In this Letter, we now analyze the generalized DM by 
coupling a ID lattice of Rydberg atoms with an optical 
cavity, where the dipole interaction between two Ryd- 
berg atoms competes with the atom-light coupling. We 
see that, while the atom-light interaction favors the SR 



phase, the atomic interaction tends to form three incom- 
pressible Rydberg solid states with filling numbers 0, 1/2, 
and 1, which destroys the polariton formation. Most im- 
portantly, we find a novel state corresponding to a "su- 
perradiant solid" (SRS) phase, where both the superradi- 
ance and crystalline orders coexist and the corresponding 
U(l) and translation symmetries are broken simultane- 
ously. Compared with the supersolid (SS) phase in an op- 
tical lattice [19-21], which breaks the same symmetries, 
the SRS is rather unique as it is induced by the non-local 
atom-light coupling and the condensation of polaritons. 
Moreover we find that, while the Solid- (0, 1) to SR phase 
transitions remain second order, both the Solid-1/2 and 
SRS to SR phase transitions become first order. 

The model - The system consists of a deep ID lattice 
with TV sites in an optical cavity. Each site is occupied 
by a single atom with the ground state |g) coupled to 
a high- lying Rydberg state |e) via a two-photon transi- 
tion process [22] . The Hamiltonian of this system in the 
rotating-wave approximation is given by 
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where u> and e are the cavity and atom transition frequen- 
cies with the detuning defined by A = lu — e. ^' is the 
single-mode creation operator of the cavity field, a, and 
bi are the boson operators representing the lower and up- 
per levels of each atom and satisfy the single-occupancy 
constraint b\bi + aja, = 1 and g s = g/\f~N is the ef- 
fective two-photon coupling. The strong dipole-dipolc 
interactions Cq between two Rydberg states is modeled 



by projectors P< 
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b\bi onto the Rydberg state, 



where only nearest interactions are considered [15]. The 
last term is the chemical potential for the total number 
of excitations A ox = iftij) + YliLi b\bi- 

The above model possesses two limiting cases. First for 
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C% = 0, the Hamiltonian (1) becomes the DM. Here we 
note that, because the ground state |g) of the atoms is not 
directly coupled to the Rydberg state |e), the so called 
no-go theorem [23] does not apply [24] , and the SR phase 
can occur. Recently, the SR phase transition has been 
observed with a superfluid atomic gas in an optical cavity 
[9] . Secondly, when g is zero this system becomes a pure 
lattice of Rydberg atoms. Then by tuning the chemical 
potential ft = /i — cj, one may derive three Rydberg solid 
states: (i) /2 < —A, it forms a Solid-0 phase with all 
the atoms staying in the ground state; (ii) —A < \x < 
2C*6 — A, half of the atoms are excited to the Rydberg 
states and form a Solid-1/2 phase because of the nearest 
neighbor repulsion; (iii) jl > 2Cq — A, all the atoms are 
excited to the Rydberg states forming a Solid-1 phase. 

Analytical Approach.~A convenient starting point is to 
consider the limit g <C Cq. In this case, one may take 
the atom-light interaction as a perturbation and use the 
strong coupling expansion (SCE) [25]. By comparing the 
energies of exciting a particle or hole with the correction 
from second order processes in the solid phase, we derive 
the following melting critical lines of the incompressible 
solid lobes [26]: 

Solid-0: p. cl = -A- 5 2 /A, 
Solid-1/2: p, c2 = -A + 5 2 /(2A), 

fi c5 = -A + 2C 6 - 5 2 /(2A - 4C 6 ), 
Solid-1: /Z c4 = -A + 5 2 /(A - 2C 6 ) + 2C 6 , 

Mc 5 - -S 2 /(A - 2C 6 ), 

which are shown as dashed lines in Fig. 2. The melt- 
ing will make the solid compressible, but it needs to be 
explored if the order in the Solid-1/2 state is actually de- 
stroyed or may be stable even in the presence of particle 
or hole excitations. 

For this purpose we introduce a variational ground 
wavcfunction [12] to gain a more physical insight into 
the phase diagram: 

\X,9) =exp(^|^)n[cos(|)6 ? t +sin(|) a i]|0>, (2) 

i 

where |0) denotes the vacuum state with all atoms in the 
ground state, and A and 9i are the variational parameters 
for the coherent cavity field and atomic fields. To find 
the ground state, we calculate the energy density £ = 
4(A,0|ir|V|A,0)/iV: 

£ = 4[— p A (sin + svclQb) + Cq cos 9 a cos 8 b 

-fiX 2 - (/} + A-C 6 )(cos0 A + cos0 s )], (3) 

where we have assumed two sublattices A and B. From 
Eq. (3) we see that while the p-term tends to enhance the 
cavity field A and the condensation of polaritons, the sec- 
ond C6-term favors the staggered order of the polaritons. 
After minimizing £ in respect to the variational parame- 
ters, we derive the phase diagram Fig. 1. The correspond- 
ing variational values are shown in TABLE I, where the 



solid phases represent the Rydberg crystals without co- 
herent cavity excitations. The usual superradiance phase 
is denoted by SR, where the atoms and light form po- 
laritons and condense with a nonzero order parameter 
(b\ai) = sin0gR,/2 and a coherent cavity field Asr. The 
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TABLE I: The variational values for different phases, where 
01 / 6> 2 in SRS phase. 



most interesting finding is that, when g ~ Cq there ex- 
ists an intermediate SRS phase between the Solid-1/2 and 
SR phases. Different from the SR phase where the polari- 
tons are excited uniformly, the SRS phase also breaks the 
translation symmetry, which shows a characteristic exci- 
tation density pA,B = (cos 6*1,2 + l)/2 that is not equal 
on the two sublattices. Therefore, both the superradi- 
ance and crystalline orders coexist in the SRS phase. In 
Fig. 1, we also see that the solid lobes shrink with de- 
creasing A. 

In bipartite lattices with nearest neighbor repulsion of 
bosons a related supcrsolid (SS) phase is known to ex- 
ist, but this phase is only stable for soft-core contact 
interaction in a ID Bosc- Hubbard (BH) model [21]. It 
is therefore surprising to find an SRS phase in our sys- 
tem which is effectively described by hard-core bosons, 
since only a single polariton can be excited in each site. 
Here we note that, different from the BH model, there is 
no direct kinetic hopping term but instead the non-local 
cavity field couples with all atoms simultaneously. This 
is the unique property that makes the appearance of the 
SRS rather nontrivial in this system. 




FIG. 1: Variational phase diagram for (a) A = 3, (b) A = 2, 
(c) A = 1, and (d) A = 0.5. All parameters are in units of 
C 6 . 

However, the above results are derived from mean-field 
calculations, so the SRS phase may not be stable when 
quantum effects are taken into account. In particular, 
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when a particle (or hole) is introduced into the Solid-1/2 
phase, the strong quantum fluctuation of the non-local 
atom-cavity coupling may drive the particle (or hole) to 
move freely on the whole lattice, which may destroy the 
staggered solid order and drive the Solid-1/2 phase di- 
rectly into the SR phase. Another possibility is that the 
SRS phase may be unstable and separates into domains 
[20] . To settle these issues and demonstrate explicitly the 
phase diagram, one needs an unbiased numerical simula- 
tion to go beyond mean-field theory. 

Quantum Monte Carlo ( QMC) simulations. -We adopt 
the high-accuracy cluster Stochastic Series Expansion al- 
gorithm [28, 29]. This algorithm is very efficient and 
has no truncation error arising from the limitation of the 
photon number. In order to distinguish different phases, 
one needs to calculate several observables: the average 
excitation density p — (N ex )/N, the structure factor 
S(q)/N = (lE^i^e^l 2 )/^ 2 , and the compressibil- 
ity k = Np{(p 2 ) - (p) 2 ). For the half filled Solid-1/2 
phase, we take the structure factor S(Q)/N with stag- 
gered order Q = tt to characterize the translational sym- 
metry breaking. 
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in a small region as shown in the enlarged region of the 
SRS phase in Fig. 2(c). In first sight it appears difficult to 
reach such narrow regions of the SRS phase experimen- 
tally. However, this changes when considering the phase 
diagram in the canonical ensemble in Fig. 2(b), where 
the SRS phase exists in a wide region Ap ps 0.1 in the 
excitation density. In Fig. 2(b) we also find that the re- 
gion of the particle-excited SRS phase is smaller than the 
hole-excited one. The phase separated (PS) regions show 
that both the Solid-1/2 and the SRS to SR phase tran- 
sitions arc first order. This differs from the mean- field 
results and may be understood by the Ginzburg Landau 
theory [27]: Because the U(l) and translation symme- 
tries are broken simultaneously, the corresponding order 
parameters S(Q)/N and (6- a.;) couple to each other and 
contribute to a sixth-order term in the free energy, which 
results in the first order phase transition. Moreover, a tri- 
critical point appears among the SRS, SR and Solid-1/2 
phases, see Fig. 2(c). 
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FIG. 3: (color online). The structure factor S(Q)/N (solid 
line) and compressibility k/10 (dash line) along the trajecto- 
ries of the lower (left panel) and upper (right panel) vertical 
cuts indicated by solid bars in Fig. 2(a). Here, the red arrows 
mark the critical points with \± = p c . For all cases = 1, 
/3 = 500, N = 100, and A = 3 (A = 2.5). The inset shows 
the corresponding average excitation density p. For simplicity, 
the small errorbars for the observables are not shown here. 



FIG. 2: (color online), (a) The phase diagram obtained by dif- 
ferent methods for Ce = 1 and A = 3 in the grand-canonical 
ensemble, (b) The phase diagram obtained by QMC in the 
canonical ensemble. Here, PS denotes the phase separation 
between the Solid-1/2, SRS and SR phase transitions, (c) The 
enlarged region of the SRS phase, where the red arrow marks 
the tricritical point among SRS, SR and Solid-1/2 phases. 



Fig. 2(a) shows the zero-temperature phase diagram 
for A = 3 and Cq = 1, where the dashed lines of SCE 
match quite well with the QMC calculations. Compared 
with the mean-field results the SRS phase is greatly sup- 
pressed, indicating that quantum fluctuations weaken the 
SRS order. Nonetheless, the SRS phase remains stable 



To provide a convincing support for the above phase di- 
agram, we now concentrate on the detailed properties of 
the phase transitions. Figure 3 shows S(Q)/N and k/10 
along the trajectories which are indicated by solid verti- 
cal lines (red) in Fig. 2(a). Here, we take T = C 6 /500, 
which is low enough to avoid the thermal effect. In both 
panels, there exist regions where a vanishing structure 
factor S(Q)/N — > is accompanied by a finite k, which 
means that the translational symmetry is not broken and 
they are compressible, which is characteristic of the SR 
phases. Regions of zero K with finite S(Q)/N correspond 
to the incompressible Solid-1/2 phases. Most signifi- 
cantly, there exist intermediate SRS phases which have 
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both finite S(Q)/N and k. The corresponding structure 
factor S(Q) /N and average density excitation p have a fi- 
nite jump and k diverges on the critical points £l = jl c . It 
signifies that the SR to SRS phase transition is first order. 
In Fig. 3, we also show the results for different A. We find 
that, while the region of the SRS phase in the left panel 
is hardly affected by A, the one in the right panel shrinks 
faster by decreasing A. This can be understood by the 
second-order atom-light interaction process which differs 
for a particle or hole excited state with corresponding 
second-order energies ~ — ( A -2C 6 ) anc ^ ~ — 7T' 

(2) (2) 

respectively. Note that \Ef'\ > \E y h '\, which means the 
particle excited SRS phase is harder to form than the hole 
excited one and more easily suppressed by decreasing A. 

We have also made a finite size scaling analysis for 
points in the SRS phase to demonstrate that the exis- 
tence of the SRS phase is robust in the thermodynamic 
limit. Figure 4 shows the corresponding scaling results in 
the SRS phase. We see that both S(Q)/N and k converge 
to finite values in the infinite size limit. 



can be excited and shared among all the atoms [18]. 

Conclusion.-We have shown that the generalized Dicke 
model of a cavity QED coupled with a Rydberg lattice 
gas displays a rich phase diagram. The competition be- 
tween the non-local atom-light coupling and the atomic 
interaction can stabilize a novel SRS phase. By imple- 
menting an unbiased QMC calculation, we find that both 
the Solid- 1/2 and SRS to SR phase transitions are first 
order, and the hole-excited SRS phase is more stable 
than the particle-excited one. This system may act as 
a new quantum simulator for the future study of quan- 
tum many-body physics. 

We acknowledge Tao Wang for helpful discus- 
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under grants No. 11704175, 10904096, NSFB un- 
der grants No. 1092009, NKBRSFC under grants No. 
2011CB921502 and by the DFG via the SFB/Transrcgio 
49. 
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FIG. 4: (color online). The finite size scaling for S(Q)/N and 
ft at Ce = 1, A = 3 and fi = 500 vs different sizes, jl — —2.905 
(—1.105) correspond to the SRS phases on the trajectories of 
the lower (upper) cuts shown in Fig. 2(a). 



Discussion. -Finally, we discuss two issues about the 
experimental implementation. First, this system could 
potentially be realized by trapping a lattice of 87 Rb 
atoms coupled to an ultrahigh finesse cavity [22]. The 
5Si/2 ground atomic state can be coupled to the 50D 5 / 2 
Rydberg excitation state [14] through a two-photon pro- 
cess via the 5J3/2 intermediate state, which can be adi- 
abatically eliminated. The resulting effective coupling 
constant <? e ff is of order 10 MHZ, which is higher than 
both the decay rates of the cavity and the Rydberg state. 
Secondly, we have assumed in the model that only a sin- 
gle Rydberg excitation is possible per site. Actually, this 
does not imply only one atom per site, but instead one 
may work with the collective degree of freedom of multi- 
ple atoms on each site, where only a single Rydberg state 
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